Competition between Electron-Phonon coupling and Spin Fluctuations in 
superconducting hole-doped BiOCuS 
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BiOCuS is a band insulator that becomes metallic upon hole doping. Superconductivity was 
recently reported in doped BiOCui-xS and attributed to spin fluctuations as a pairing mechanism. 
Based on first principles calculations of the electron-phonon coupling, we argue that the latter is very 
strong in this material, and probably drives superconductivity, which is however strongly depressed 
by the proximity to magnetism. We find however that BiOCui-^cS is a quite unique compound where 
both a conventional phonon-driven and an unconventional triplet superconductivity are possible, 
and compete with each other. We argue that, in this material, it should be possible to switch from 
conventional to unconventional superconductivity by varying such parameters as doping or pressure. 

PACS numbers: 63.20.Kd,74.20.Pq,74.20.Mn,74.70.Xa 



The study of spin fluctuations as superconducting me- 
diators dates back to the sixties jii^ however, in contrast 
to the electron-phonon (EP) interaction, for which a de- 
tailed first-principles theory has been developed in the 
last twenty years, a quantitative theory is still lacking. In 
several materials where at some point ferromagnetic spin 
fluctuations (paramagnons) were considered as potential 
pairing agents, such as ZrZn2,'^ MgCNia,^ or Pd metal, ^ 
phonon and spin fluctuations contributions either cancel, 
rendering the material non-superconducting (ZrZn2, Pd), 
or the latter substantially decreases the superconducting 
transition temperature. 

Recently, superconductivity with rc=5.8 K has been 
discovered in hole-doped BiOCuo.gS.^ BiOCuS crystal- 
lizes in the ZrCuAsSi-type structure, isostructural to the 
1111-family of Fe-based superconductors, with Cu-S lay- 
ers playing the role of Fe-As layers. While Cu-S hy- 
bridized dp bands per se are rather similar to the Fe- 
As bands in Fe-pnictides, the different electronic filling 
brings about very different properties in the two systems. 
The stoichiometric BiOCuS is in fact a band insulator 
with the Cu being in the d^*^ electronic configuration;^^- 
With hole-doping it displays both a strong tendency to 
itinerant (ferro)magnetism, and a spectacularly strong 
EP coupling, hinting to unconventional, triplet p-wave^ 
and conventional, singlet s-wave superconductivity, re- 
spectively. 

In this paper, we study the interplay between these 
competing instabilities, using first-principles calculations 
of BiOCui_2;S as a function of doping and Stoner param- 
eter, which we use as a proxy for the tendency to mag- 
netism. We find that, as the EP coupling is spectacularly 
strong, it is likely that a conventional superconductiv- 
ity, even though depressed by spin fiuctuations, is more 
stable than an unconventional {e.g. p— wave) one. It ap- 
pears though that a small variation of parameters can re- 
verse the situation and bring triplet superconductivity or 
long-range magnetism. We identify two large regions in 
the parameter space where, respectively, ferromagnetism 



(FM) or conventional s— wave superconductivity are the 
ground states, with an intermediate region where no FM 
long range order is predicted, yet spin fluctuations are 
strong enough to destroy the s-wave superconductivity 
and possibly stabilize a triplet state. 

We perform calculations in the linear-response approx- 
imation for the EP interaction, and in the local spin den- 
sity functional version of the random-phase approxima- 
tion (RPA) for spin fiuctuations, as described below; dop- 
ing is treated in the rigid band approximation (RBA)i^ 

The generalized-gradient approximation (GGA) band 
structure and partial electronic density of states (DOS) 
are shown in Fig. [TJ in agreement with previous calcu- 
lations ,^i^i^^ we find that the stoichiometric compound 
is a semiconductori^ with an indirect gap of A « 0.5 
eV (GGA); the top of the valence band occurs along the 
T — M line, and we choose it as the zero of the energy 
in the following. The electronic structure in an energy 
range ~ 7 eV below the top of the valence band in BiO- 
CuS is derived from Cu d and S p states (see top panel of 
Fig. [1]). The Cu d states are centered around ^ —2 eV. 
They hybridize strongly with the S p states, forming anti- 
bonding bands within ^ 1 eV below the semiconducting 
gap. The EP matrix element is large for these bands, as 
the electronic states are very sensitive to ionic displace- 
ments. On the contrary, the deeper, non-bonding, Cu 
d bands, centered around ^ —3 eV, are less sensitive to 
the Cu-S hopping parameters and exhibit a much weaker 
EP interaction. The tendency to magnetism is instead 
strong throughout the entire Cu d band, since the Stoner 
parameter of Cu is large {Icu ~0.9 eV). 

In pure BiOCuS, Cu is in a nominal d^*^ state and thus 
not magnetic. Doping with holes, for x < 0.5, shifts the 
Fermi level into a doubly-degenerate band, with domi- 
nant Cu dxz, dyz and S Px, Py characters and large DOS. 
The bottom panel of Fig. [T] shows a blow-up of the band 
structure in the energy range relevant for superconduc- 
tivity. The dotted and dashed lines indicate respectively 
the position of the Fermi level for x = 0.1, corresponding 



to the doping for which superconductivity was observed 
in Ref. [B, and x — 0.5, which is the highest doping con- 
sidered in our RBA study. 

If we could shift the Fermi level further down, so as to 
cut the band structure at ^ —1.4 eV (dash-dotted line 
in the top panel of Fig. [1} , we would find a striking sim- 
ilarity with the familiar low-energy electronic structure 
of Fe-pnictides, with the xz,yz hole and electron pock- 
ets, centered at F and M respectively. The DOS and 
the p — d hybridization here are small, thus the tendency 
to antiferro- (rather than ferro-) magnetism, and low EP 
interaction. This is indeed what first-principles calcula- 
tions find in Fe-pnictidesJ^ 

We now go back to discuss the behavior of BiOCui_a;S 
for X < 0.5, using the bottom panel of Fig. [T] For x > 0.1, 
we find the ground state of the system is FM, both in 
the local spin density approximation (LSDA) and in the 
GGA.— This can be understood in terms of the Stoner 
criterion for FM, INq > 1. When holes are introduced 
into the system, and the Cu charge state is being reduced 
from to cP , the Fermi level moves into a flat region of 
the band structure, which gives rise to a high peak in the 
DOS (Nq = 2.1 st/eV spin). Since the Stoner parameter 
of atomic Cu is Icu ~ 0.9 eV, the latter value of the 
DOS is well above the Stoner criterion for FM, INq > 1. 
In BiOCui_a;S the actual value of / < Icu, due to the 
Cu-S hybridization. It can be estimated from the split- 
ting AE = ml between majority and minority bands in 
the FM state, where m is the value of the self-consistent 
magnetic moment. We find m < 0.1 for all dopings con- 
sidered, and / = 0.53 eV in LSDA and / = 0.67 eV 
in GGA, independent of doping. So far, however, experi- 
ments have seen no trace of static magnetism; this is con- 
sistent with the tendency of LSDA calculations to overes- 
timate the tendency to itinerant magnetism with respect 
to experiment in the vicinity of a magnetic quantum crit- 
ical point (QCP), where the system exhibits strong spin 
fiuctuations 1^ We will return to this issue in more detail, 
after discussing the results for the EP interaction. 

Fig. [2] summarizes the EP properties of the hole- 
doped BiOCui_a;S. The partial phonon density of states 
(PDOS) of the undoped compound extends up to 65 
meV; vibrations of the Bi-0 layers are concentrated at 
energies < 20 meV, while modes involving the Cu-S lay- 
ers are found at higher energies. The S atoms give rise to 
a very broad feature in the PDOS, from 40 to 65 meV. 
Using this phonon spectrum, we calculate the Eliashberg 
spectral function a'^F{uj) of the hole-doped BiOCui_a;S: 



^ ^ No ^ 




R X M 



ro 2 4 



FIG. 1: (color online), top Band structure of BiOCuS, shaded 
according to the partial Cu dxz+yz (left) and S Px+y (right) 
characters: the continuous and dashed-dotted lines mark re- 
spectively the position of the Fermi level in the undoped com- 
pound and that corresponding to the filling of Fe-pnictides 
(see text); the corresponding DOS is also shown, bottom: a 
blow-up of the low-energy band structure; the dashed and 
dotted line mark the position of the Fermi level, correspond- 
ing to a hole doping a; = 0.1 and x — 0.5, in RBA. 
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FIG. 2: (color online). From top to bottom: Partial Phonon 
density of States (PDOS), Eliashberg spectral function for 
X — 0.1, in RBA, and (inset) ratio between the coupling con- 
stant and the DOS as a function of doping. 
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evaluating the average of the EP matrix elements 
on the Fermi surface S{e^), obtained by a rigid-band shift 
corresponding to the doping level. From the Eliashberg 
function we calculate the EP coupling constant: 

= 2 dna'^F{n)/n. 

For all dopings x < 0.5, we find that only two groups of 



phonon modes, corresponding to the out-of-plane vibra- 
tions of the Cu-S layers, have sizable EP matrix elements 
■9k k+q- these give rise to two narrow peaks in a^F^uj), 
centered at 33 meV and 50 meV. The lower panel of Fig. [2] 
shows an example of a^F^oj) for x = 0.1. 

Since the shape of the Eliashberg function does not 
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depend on x for all dopings considered, the total EP 
coupling depends on doping only through the value of 
the density of states at the Fermi level, iVo. We thus 
rewrite Acp as Acp — NoVcp- As the inset of Fig.[5]shows, 
Vcp ~ 0.9 eV spin f.u. at all dopings for a; < 0.5. For com- 
parison, Vcp = 0.1 cV spin f.u. in LaOFeAs and Vcp — 0.3 
eV spin f.u. in Pd (i.e. in metals where the lattice plays 
a minor role compared to spin fluctuations) while it is 
much larger in good EP superconductors: Vcp = 2.5 eV 
spin f.u. in MgB2 or Vcp = 6.6 eV spin f.u. in Pb. 

For X = 0.1, A^o = 1-93 st/eV spin f.u., Acp = 1.74 and 
the logarithmically averaged phonon frequency wiog = 
263 K. This EP interaction would then give rise to a 
Tc of 33 K, assuming a typical value for the Coulomb 
pseudopotential, /i* = 0.1. 

This is much larger than the experimental value Tc = 
5.8 K,^ which would correspond to Acp = 0.6.-^- A factor 
three discrepancy is well above the typical uncertainty 
of Tc in similar calculations, stemming from the compu- 
tational uncertainty on Acp, typically 10%, or from the 
uncertainty of ^* . 

The most straightforward explanation, in the present 
case, is a suppression of phonon-mediated pairing by 
strong paramagnons, due to proximity to a FM QCP. We 
now estimate this effect, using the RPA. Let A^f be the 
coupling to spin fluctuations in the singlet channel; the 
effect of paramagnons is to suppress superconductivity in 
the singlet channel by depressing the effective coupling 
constant (Aa = Acp — A^f) and increasing the effective 
mass of the carriers by the factor 1 -|- A^ = 1 -I- Acp -I- A|f. 
This effect has been studied in Ref.[l^ where the following 
expression for Tc was derived (and verified by comparison 
with numerical solutions of the Eliashberg equations) : 



Tc = 



1.45 



exp • 



-(1 + Az) 



Aa-m*(1 + 0.5^) 



(1) 



Here we assume for simplicity that the characteristic fre- 
quencies of phonons and electrons are the same. Eq. [1] 
can also be generalized to triplet superconductivity, with 
the substitution: Aa A*f ; Xz ^ Xz — ^cp + ''^sf' where 
A*f = iAgf is the coupling to spin fluctuations in the 
triplet channel.— Eq. [T] gives an appreciable Tc only if 
the denominator in the exponential is positive. For small 
fj,*, this is the case, when Aa > 0. We therefore use Aa to 
define the phase diagram of hole-doped BiOCui_a;S: us- 
ing the RBA, we take Acp (a;) = VcpNo{x), where No{x) is 
the DOS at the Fermi level at doping x. For the coupling 
to spin fluctuations we use the following expression: 
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where / is, in the LDA parlance, the Stoner parameter i2i 
Eq. [2]is similar to the well-known expression for the spin 
fluctuations induced interaction in the singlet channel,— 
averaged over the Fermi surface. Note that in the triplet 
channel the spin fluctuations interaction is three times 



smaller, and also the averaging for both Agf and Acp is 
performed with a weighting factor WF(k) • vp^li.'). 

A well-known LDA problem is that, due to its mean- 
field character, it overestimates the tendency to static 
magnetism. This can be corrected by introducing a phe- 
nomenological Stoner /, reduced from its LDA value. In 
this spirit, in the following, we treat / as a free parame- 
ter, and plot the phase diagram of BiOCuS in the {x, I) 
space. The results are shown in Fig. [3] 



2.0 




X 



A 

1.0 
0.0 



_l I I I I I I J 

0.10 0.20 0.30 0.40 0.50 
X 



FIG. 3: (color online). Phase diagram of BiOCui-^S, de- 
fined by Aa as a function of doping (x) and Stoner param- 
eter I. whose value is represented by the color scale. The 
two horizontal dashed lines correspond to /lda=0.53 eV and 
-^GGA=0.67 eV. The vertical dashed line indicates the doping 
for which superconductivity was observed in Ref. @. In the 
region (FM) the system shows a FM instability, defined by 
the condition (A^o^ > 1); elsewhere the system is paramag- 
netic (FM). Below the bold line (which marks the condition 
Tc = Tc) the ground state is a conventional singlet supercon- 
ductor. Above the bold line a triplet superconducting state 
is more stable. The isolines Aa = 0.6 and A^f — 0.6 indicate 
the values of I, x which reproduce the experimental Tc = 5.8 
K of Ref. ^ in the singlet and triplet channel respectively. 

If Aa 3> M* a conventional EP superconductivity, al- 
beit depressed by spin fluctuations, is a stable zero- 
temperature ground state. As the Stoner parameter is 
increased Aa goes down, and a competing instability 
against a triplet state emerges when the critical temper- 
ature in the singlet channel T^, defined by Eq. [U be- 
comes equal to that in the triplet channel (T*). Finally, 
as the tendency to magnetism is increased even further, 
the Stoner criterion A^o^ > 1 is satisfied, and the system 
becomes ferromagnetic (Fig. [3]). 

One can see that, had we used the LDA or GGA value 
for /, for dopings close to x = 0.1, we would have found 
BiOCuS inside the FM region. However, at x = 0.1 ex- 
periments see no trace of static FM order, a sign of in- 
adequacy of the mean field character of magnetism in 
LSDA. Reducing the LDA value of / to Jeff = 0.51 eV 
suppresses the magnetic instability at x = 0.1; a reduc- 
tion to Jeff = 0.39 eV brings the estimated triplet Tc into 
agreement with the experimental one, and a reduction 



4 



to 0.25 eV does the same with the conventional singlet 
Tc- For typical itinerant magnets renormalizing Ilda by 
'-^ 30 — 40% provides reasonable agreement with the ex- 
perimental magnetic susceptibilities,^^ in the same ball- 
park as the reduction introduced above. 

In other words, BiOCuo.gS is a unique example where 
a spin fluctuations driven triplet superconductivity is 
nearly degenerate with the phonon-driven singlet super- 
conductivity, and the critical temperature is sizable for 
both symmetries. Given that the actual BiOCuo.gS sam- 
ples are rather dirty, one may conjecture that samples 
studied in Ref. d are on the conventional side of the phase 
diagram, but the fact that superconductivity appears to 
be so difficult to reproduce may be due to the fact that 
slightly different samples may appear outside of the sta- 
bility range of singlet pairing in the phase diagram in 
Fig. [3l In principle, one can use pressure and doping, 
which control / and A'o respectively, to move around in- 
tentionally in the proposed phase diagram. 

This tunability comes about because of the combina- 
tion of two factors: an exceptionally strong EP inter- 



action in the singlet channel that is essentially canceled 
out in the triplet channel, and a strong spin fluctuations 
coupling that competes with EP interaction in the sin- 
glet channel. The occurrence of these two large coupling 
constants can be seen as the result of three concurring el- 
ements: a strong d—p hybridization, that causes large EP 
matrix elements; the large value of the Stoner parameter 
of Cu, that causes a strong tendency to magnetism; and, 
finally, the presence of a large peak in the electronic DOS, 
which favors FM and enhances the coupling constants for 
superconductivity both in the singlet and triplet channel. 
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